Cascades of energy and helicity in the GOY shell model of turbulence 
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The effect of extreme hyperviscous damping, uk^,p = oo is studied numerically in the GOY 
shell model of turbulence. It has resently been demonstrated [Leveque and She, Phys. Rev. Lett, 
75,2690 (1995)] that the inertial range scaling in the GOY model is non-universal and depending on 
the viscous damping. The present study shows that the deviation from Kolmogorov scaling is due 
to the cascade of the second inviscid invariant. This invariant is non-positive definite and in this 
sense analogous to the helicity of 3D turbulent flow. 
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The invicid invariants, like the energy, of the Navier-Stokes equation are important quantities determining the 
dynamics of turbulent flow. The major difference between 2D and 3D flow originates from existence of a second 
inviscid invariant, enstrophy, in 2D - absent in 3D flow. This gives strong constraints on the 2D flow dynamics. 
Conservation of enstrophy leads to forward cascade of enstrophy and backward cascade of energy in 2D, giving rise 
to large-scale coherent structures in the flow. In 3D there also exists a second quadratic inviscid invariant, namely 
the helicity, defined as the integral of the scalar product of the velocity and the vorticity. It has been proposed to 
be important for 3D turbulence 0]. It was shown by Kraichnan 0] that the interaction between waves of opposite 
helicity is stronger than interactions between waves of same helicity. Consequently, helical flow might slow down the 
Q> ' cascade of energy to the dissipation range . 

Since the spectral helicity density grows with a factor fc, the modulus of the wave- vector, over the the energy density 
the helicity dissipation would, in the R — ► oo limit, prevent energy dissipation. The situation would resemble that 
of 2D turbulence where cascade of enstrophy prevents forward cascade of energy. In 2D turbulence the conservation 
of both J dku(k) 2 and J dkk 2 u(k) 2 implies that forward cascade of enstrophy must be accompanied by a backward 
0^ . cascade of energy. This strong constraint does, however, not apply in the case of helicity transfer since helicity is not 
positive definite. So we only have the weaker assessment mentioned above that the efficiency of energy dissipation 
could be depending on the non-linear helicity transfer and the helicity dissipation. 

Considerable interest has lately been given to the GOY model of turbulence, introduced by Gletzer and examined 
by Yamada and Ohkitani jjj. Comprehensive lists of resent references on the GOY model can be found in refs. H|). 
' The behavior of the helicity in different shell models have resently been investigated ||J^] . 

The GOY model is a simplified reduced wave-number analog to the spectral Navier-Stokes equation. The spectral 
domain is represented as shells, each of which is defined by a wavenumber k n = ko\ n , where A is a scaling parameter 
. defining the shell spacing. There are 2N degrees of freedom, where iV is the number of shells, namely the generalized 
' complex shell velocities, u n for n = 1,N. The dynamical equation for the shell velocities is, 

H [ be 

u n = ik n {au* n+2 u* n+1 + -r< + i<_i + T2<-i<_ 2 ) - vk v n u n + /<5„ ; „ , (1) 

where the first term represents the non-linear wave interaction or advection, the second term is the dissipation, and 
the third term the forcing, where no is some small wavenumber. Throughout this paper the standard 3D GOY model 
parameter values, A = 2, ko = A -4 , a = 1, b = c = — 1/2 is used. The GOY model in this form contains no information 
about phases between waves, thus there cannot be assigned a flow field in real space. 

The model has two conserved integrals, in the case of no forcing and no dissipation [y = / = 0) (inviscid invariants). 
These are, E — 5X^=1 \ u n\ 2 and H — A l) n k n \un\ 2 which corresponds to the conservation of energy and 

a second non-positive definite quantity interpreted as analogous to helicity Q for the Navier-Stokes equation of 3D 
turbulence, hereafter referred to as the helicity. It should be stressed that this analogy is only in the sense that both 
quantities are non-positive definite. In this model each shell is maximally helical with alternating sign, since the 
numerical value of the helicity density is k times the energy density. 

The forcing in (|l|) is applied at a small wave number and the dissipation dominates at large wave-numbers, so 
we can define an inertial range where the non-linear energy cascading terms dominate and the Kolmogorov scaling 
arguments apply. For an energy cascade we have the "Kolmogorov scaling" of the shell- velocities, \u\ ~ 77 1 / 3 fc" 1 / 3 , 
and for a helicity cascade we have the "helicity scaling", \u\ ~ rf-^k^ 2 ^, where 77 and 77 are the mean dissipation 
per unit time of energy and helicity respectively. The model has both the Kolmogorov scaling, u n = k n g(n) and 
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the helicity scaling, u n = (— l) n k n 2 ^ 3 g(n) as unstable fixed points in the unforced and inviscid case. The function, 
g(n)=g(n+3), is any mod(3) function. The mod(3) symmetry is an artifact of the GOY model discussed in detail 
in ref. ||. It will become important in the following. The Kolmogorov fixed point plays an important role for the 
behavior of the GOY model, with forcing and dissipation, in the sence that the phase space trajectory of the shell 
velocities seems to "curl around" this point with the average values of the velocities close to the fixed point values. 
It was shown in a numerical study by Leveque and She || that the inertial range scaling in the GOY model is not 
universal and depending on the form of the viscous damping. They attributed this non-universality to the reflection of 
energy flow from the viscous subrange back into the inertial range in the case of hyperviscosity (p > 2). By studying 
the extreme case p = oo, I suggest a different explanation for the non-universality of the inertial range scaling, namely 
that cascade of helicity blocks the cascade of energy and thus changes the scaling. This effect is probably specific to 
the GOY model, and a consequence of the specific odd-even asymmetry of the helicity in the GOY model. 

The numerical study is performed on 2 versions of the GOY model. Firstly, the usual model with normal (p = 2) 
dissipation and secondly, with viscosity only applied on outermost shell, corresponding to p = oo hyper- viscosity. The 
forcing in (|l|) is in both cases taken to be / = fo/u* lQ where fg is a constant. This gives a constant input of energy 
(and helicity) per unit time. 

Figure 1 shows the result of a numerical calculation for fo = (1 + i) x 10~ 3 ,rto = 4 and v = 1CP 5 . The model 
has 20 shells and is run for 3000 time units after a 1000 time units spinup starting from the (unstable) Kolmogorov 
fixed point. First panel shows log A ((|w„|)) as a function of the shell number n — log\(k n ), where (.) denotes temporal 
average. The spectral slope is close to the Kolmogorov scaling, shown by the line, thus the model shows an energy 
cascade. 
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FIG. 1. The standard GOY model with / = (1 + i) X 1CT 3 , no = 4 and v = 10" 5 . The model has 20 shells and is as before 
run for 3000 time units after a 1000 time units spinup starting from the (unstable) Kolmogorov fixed point. log A ((|u,i|)) is 
shown as a function of the shell number n = log\(k n ). The line indicates the Kolmogorov scaling with spectral slope of —1/3. 



The sources and sinks for energy and helicity are 

E = - vklul) (2) 

n 

and 

6 = E((- 1 )" fc ™<^ ~ (3) 

a 

This means that helicity of opposite signs is dissipated at every second shell and therefore total helicity is produced 
for odd numbered shells in the dissipation range. 

The inertial range flow is dominated by the non-linear transfer of the conserved quantities. The key assumption 
by Kolmogorov (1941) is that the inertial range flow does not depend on the specific form of the small scale viscous 
dissipation. This is not the case for the GOY model. In order illustrate this point a study of a slightly modified 
GOY model is presented here. The dissipation is taken to be active only on the outermost shell, number N = 20, 
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this corresponds to hyperviscosity in the limit, p =4> oo, with v = v§k~^f and Vq held constant. With this choice of 
dissipation the behavior of the GOY model changes dramatically even though the non-linear terms in the governing 
equation remains unchanged. 

Figure 2 shows the result of a numerical integration with the modified GOY model where /o = (1 + i) x 10~ 8 
and Vq = 10~ 7 fc| = 429.5. First panel shows the velocity spectrum where the mod(3) symmetry of the shell model 
becomes dominant. The scaling behavior of shells numbered 10,13,16,19 is different from the rest. The line in the 
figure has a slope of —2/3 corresponding to the cascade of helicity rather than energy. 

However, if U n = \u(n — l)u(n)u(n + l)! 1 ^ 3 is considered the mod(3) symmetry is eliminated and the inertial 
range scaling approximately reemerges as an arithmetic mean of the scaling behaviors for the shells 3n, 3n + 1, 3n + 2 
(figure 2, second panel). The line in the figure has a slope of —1/2 corresponding to equipartitioning of helicity. 
The hyper-viscosity of the model, only pulling out helicity through the (18,19,20) triade, cannot maintain a helicity 
cascade. 
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FIG. 2. The modified GOY model with f = (1+ i) X 10" 8 and v = 10~ 7 fc|, = 429.5. This model has likewise 20 shells 
and is run for 3000 time units. First panel is the same as in Figure 1. The line in the figure has a slope of —2/3 corresponding 
to the cascade of helicity. Second panel shows the spectrum of U„ = \u(n — l)u(n)u(n + l)] 1 ^ 3 . The line has a slope of —1/2 
corresponding to equipartition of helicity. 

The result can be interpreted as a spectral bump at the end of the spectrum with an approximate equipartition of 
helicity. Similar results have been found by Borue and Orszag || for numerical simulations of the 3D Navier-Stokes 
equation. In this case the accumulation seems to be an effect of energy not being able to be transferred across the 
ultra-violet cutoff. 

In the case (/ = v — 0) the model has, discarding the boundary effects, besides the two scaling fixed points a 
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periodic solution. It is easy to verify that the following satisfies the dynamical equation (1): 

U3„-2(t) = sk^X- 3 " 

«3n(t) = A 3 "T^e^- aia2t 

(4) 

with 

a 1 = s\- A (l + b\ 2 + c\- 3 -> +1 ) 
a 2 = sA 37 ~ 3 (l + bX- 3 "'- 1 + cA- 37+1 ), 

(5) 

and s is an arbitrary constant. The scaling parameter, 7, is related to the scaling fixed point by 7 = —(a + l)/3, 
where z = X a is a solution to 1 + 6z + cz 2 = 0, thus a generator of one of the conserved quantities ITS). So there are 
two values possible for 7; 7 = — (2 + wr/ log(A))/3 corresponding to the generator, z — —2, of helicity, or the fluxless 
fixed point of the GOY model, and 7 = —1/3 corresponding to the generator, z = 1, of energy, or the Kolmogorov 
fixed point. 

There is in this periodic solution a complete phase locking of all the shells. The energy and helicity fluxes are both 
zero in the periodic solution, so in some sence it corresponds to the fluxless fixed point of the GOY model. 

The shell-velocities of shells 8, 11, 14, 17, 20 arc out of phase with those of shells 9, 12, 15, 18, while shell- velocities 
of shells 10, 13, 16, 19 are almost constant. From this it seems as if the periodic solution plays the same role for this 
modified GOY model as does the Kolmogorov fixed point for the GOY model. This could indicate the existence of 
an (unstable) limit cycle in the modified GOY model. 

The dissipation of helicity is irregular, associated with the bursts in shell number 20 (and all other shells), the 
energy dissipation is completely blocked, resulting in a steady increase of energy concentrated on shells 10,13,16,19. 
The mechanism for preventing the forward cascade of energy in this model is analogous to the 2D case, where the 
forward cascade of enstrophy governs the dynamics. 

The system shows no signs of approaching a statistical equilibrium simply because the helicity can grow to arbitrarily 
large negative values (dominated by shell number 19). The reason for this can be understood by examining the 
dissipative energy balance. In the case of normal viscosity it follows from (0) that the phase space trajectory will be 
attracted to the hyper-ellipsoid given by ^ vk^\u n \ 2 = f which is a compact 2N — 1 dimensional object. In the p = 00 
case this object will be the 2N — 1 dimensional "hyper-cylinder" defined by |«2o| — //^lo)- This is not compact 
and ergodicity does not apply, thus the p — > 00 is a singular limit and no statistical equilibrium can be reached. The 
relative change in energy over the integration is of the order 10 -2 such that the system is in a quasi-equilibrium state 
where statistical equilibrium is reestablished in the limit where the energy injection rate goes to zero. Taking an even 
smaller forcing does not change the statistics, thus the statistical timeaveraging is meaningful. 

If the energy cascade is effectively blocked the energy should be pulled out of the system by a drag at the small 
wave numbers corresponding to to backward cascade in 2D. This is not seen, so an invers energy cascade cannot be 
established. The model was also run with the usual dissipation but only active on every second shell. The result of 
this run was essentially the same as for the modified GOY model. 

In conclusion we see that the non-linear transfer in the GOY model depends crucially on the dissipation properties 
for both conserved quantities, energy and helicity. It is still an open question how much the GOY model reflects the 
dynamics of the Navier-Stokes equation. These findings certainly indicates that the model is to restricted in some 
sence to represent real flow. The dissipative term introduced in the modified model is certainly not very realistic, in 
real flow helicity of both signs will be dissipated, which it is not in the modified model presented here. The findings 
from this simple model might indicate that in the inertial range flow the conservation of helicity could block the 
energy cascade and thus alter the Kolmogorov fc -5 / 3 scaling. 
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